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Abstract 

An algorithm is constructed to derive a small momentum expansion for two-loop 
two-point diagrams in all cases where, due to the presence of physical thresholds, 
there are singularities at zero external momentum. The coefficients of this "zero- 
threshold" expansion are calculated analytically for arbitrary masses. Numerical 
examples, using diagrams occurring in the Standard Model, illustrate the conver- 
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1. Introduction 



The increasing precision of experiments |] testing the Standard Model (SM) of elec- 
troweak interactions will demand a matching precision from the theoretical predictions. 
The present LEP/SLC experiments measure the Z mass, its total and partial widths, 
and various asymmetries with high accuracy. Theoretical predictions for these quantities 
depend on the unknown top quark and Higgs boson masses so that some freedom is left to 
fit theory to experiment. When knowledge of the top mass from direct detection becomes 
available the possibility to adjust in the theoretical predictions will be narrowed and 
the tests will become more stringent. At present there are hints of a value for the top 
mass of 174 ±1011^ GeV 0. 

A similar situation holds for the W mass. At present its measured and predicted 
(from the Fermi constant Gp) values are in agreement when suitable vnt and Mu values 
are chosen. However the expected higher accuracy on My^ from LEP2 measurements 
will test the theory more severely. 

When at some point in the future discrepancies arise between theoretical predictions 
and experiment it will be tempting to invoke extensions of the SM. This can only be done 
convincingly when SM predictions really fail. This brings us to the question how accurate 
the present theoretical calculations are. 

The one-loop predictions of the SM have been studied by many people and are incor- 
porated in various fitting programs in use at LEPl. Some special second order effects are 
often included as well. Although in principle this should be well under control, in practice 
there are subtle differences in implementing the theory, which may cause small deviations 
in theoretical predictions. An analysis of this situation is at the moment in progress [Q. 

For some quantities one-loop calculations are definitely not adequate and thus the 
dominant parts of some second order contributions should be added. The clearest ex- 
amples are various QED radiative corrections to the Z-line shape and to small angle 
Bhabha scattering 0? 0- The former are known to the desired accuracy, the latter 
not. For small angle Bhabha scattering one could for instance be interested in double 
box diagrams |jlO|, which involve massless and massive particles. Examples of second 
order corrections in the non-QED part of the theory are contributions to the p parameter 
III , |12] and to Z ^ bb decay, for a review see |]l3l . Here one encounters two- loop diagrams 



with many massive particles. The present calculations consider some limiting situations 
of heavy top mass or heavy Higgs mass, which allow the neglect of Mz or Mw masses. 

In the future these types of approximations will become less acceptable. Therefore it 
seems unavoidable that physics issues will dictate the study of exact two-loop corrections 
to the electroweak theory. 

This is why during the last few years a number of studies have been made into the 
question of massive two-loop diagrams, the purely massless case being much better un- 
derstood. It is a sensible strategy to study firstly the simplest two-loop diagrams, the 
self-energy diagrams. Not surprisingly, most papers deal with those diagrams. From a 
physical point of view this can also be justified, since at the one-loop level these diagrams 
are usually the most relevant ones. Since it has been shown that all two-loop self-energy 
diagrams can be reduced to scalar diagrams ||T^ (see also [|T^), it is sufficient to consider 
only the latter ones. 
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The general case of massive two-loop diagrams occurs to be very complicated, even 
for two-point (self-energy) functions. Exact results in terms of known special functions 
are known for a restricted number of special cases only (some of them can be found 
in |16|, |l^ |l8l). In all these examples at least one of the masses is equal to zero. A 
difficult problem, unknown at the one-loop level, arises when a diagram contains a three- 
particle cut where all particles are massive. There are some arguments that this result 
expanded around four dimensions is for general fc^ not expressible in terms of known 



functions like polylogarithms [|T^ (see also in [16| where a special case of such a diagram 
was investigated). 

To be able to deal with diagrams that have not been evaluated exactly, a number 
of approaches have been developed, and some useful results have been obtained. One 
method is a reduction to a two-dimensional integral representation, both for convergent 
20 1 and divergent ||21| case^. A number of one-dimensional integral representations with 
more complicated integrands are available [^]. Some other methods [^, ^ can also 



be employed, although they involve representations with a larger number of parametric 
integrals. In all these cases, the final result is obtained by numerical integration. 

When one is exclusively interested in numbers and less in the analytic structure the 
above method is sufficient, although there is always the question whether the numerical 
program works correctly in all cases. When one wants to understand the analytic structure 
and also wants to have another, independent method, series expansions in the variables of 
the problem, i.e. k"^ and the masses, provide a good alternative. Moreover, this approach 
can be very useful when some of the masses are unknown, because it gives analytic 
expressions for the coefficients of the expansion. 

For a general massive scalar two-loop self-energy diagram the following expansions in 
A;^ have been studied. For small k"^ an algorithm was developed for the construction of the 
Taylor series coefficients (the three-point case was considered in [^0) . The expansion 
is valid below the smallest physical threshold. For large k"^ again an algorithm for the 
evaluation of the expansion coefficients was developed . This is based on general 
results on asymptotic expansions of Feynman diagrams where the coefficients are 
conveniently expressed in terms of more simple diagrams of a certain type. This expansion 
is valid above the largest physical threshold. Although general explicit expressions for 
coefficients of arbitrary order have not been obtained in such a way, for some special 
diagrams multiple series in k'^ and the masses, with closed expressions for the coefficients, 
are now available These multiple series expansions converge in certain regions and 



correspond to generalized hypergeometric functions. 

From this discussion it is clear that there are k^ regions, where no expansion is avail- 
able. Namely, the behaviour above the smallest but below the largest physical threshold 
needs to be described. It is the purpose of the present paper to start to fill this gap. 
Clearly filling the most general gap would be the best. But we are going to start by filling 
a smaller gap, corresponding to the cases when the lowest physical threshold is zero, so 
besides massive particles some massless particles are present in the diagram. These are 
cases when an ordinary Taylor expansion does not work because of infrared singular- 

^ Analogous integral representations can be written also for three-point two- loop diagrams p^ . 
■^In ref. p7|] , conformal mapping and Pade approximations were also employed for numerical study of 
the behaviour above the threshold. 
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ities connected with the hmit 0. By doing this, we cover all remaining holes in the 

small momentum expansion. Moreover, we expect that this specific case may be helpful 
for the most general problem of describing threshold behaviour. 

In this paper, we are going to examine the region between the k"^ = threshold and the 
first non-zero threshold. It will be shown that general results on asymptotic expansions 
of Feynman diagrams can also be employed in this case. 

The actual outline of the paper is as follows. In section 2 the algorithm for the 
expansion is constructed from the general theorem. All possible zero-threshold cases are 
considered. The actual calculation of several coefficients is presented in section 3. Also 
a comparison with some known analytical results is made. Section 4 contains numerical 
comparisons for some complicated, analytically unknown cases. Conclusions are given in 
section 5. In the appendices we present some formulae for two-loop vacuum integrals and 
discuss the evaluation of massless integrals with numerators. 



2. Constructing the expansion 



As in the previous papers |l26| , we shall consider here two-loop Feynman integrals 



corresponding to the self-energy diagram of Fig. la, denoting them as 

d> d"g 



■J i{^^h {mih k) = I I n-2 n-s n^4 n^s > (1) 



where n = 4 — 2^ is the space-time dimension (in the framework of dimensional regular- 
ization |3^), and are the powers of the denominators Di = — m'^ + iO, pi being the 
momentum of the corresponding line [pi are constructed from the loop momenta p and q 
and the external momentum k). We are mainly interested in all z/'s being integer. The 
cases when some of the z/'s are zero correspond to reducing lines in Fig. la to points. In 
such a way, self-energy diagrams with four or three internal lines can also be described. 
Moreover, by trivial decomposition (partial fractioning) of the first and the fourth denom- 
inators (provided that z/i and z/4 are integer) we can reduce the integral corresponding to 
Fig. lb to the integrals ([^) with ui or 1/4 equal to zero (such a decomposition is required 
only if mi 7^ m^). So, in the general case of self-energy diagrams (with integer z/'s) it is 
sufficient to consider only the integrals (|I|). 

In general, the diagram in Fig. la has four physical thresholds. Two of them corre- 
spond to two-particle cuts, 

fc^ = (mi + 7714)^ and k^ = (m2 + m^Y, (2) 
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while the other two are connected with three-particle cuts, 



k'^ = {nil + ^3 + TTT'b)^ 



and 



k"^ = (m2 + + m^) 



(3) 



In special cases, some of these threshold values can be equal. 

If all the thresholds (0)-® are not equal to zero, the integral (0) is an analytic function 
in a vicinity of fc^ = 0, so that the asymptotic expansion is given by its Taylor series in 
k"^. It is convergent when k"^ is less than the value of the first physical threshold. The 
general case of such an expansion of two-point functions was considered in pH], and some 



results for three-point functions can be found in [27|. The situation is different, however, 
when some of the physical thresholds (^-(^) become equal to zero. In this case, there 
can be non-polynomial terms in fc^ as fc^ — > 0. Below we shall see that they appear as 
ln(— fc^) terms as ^ 4. These are the "zero-threshold" cases that will be considered in 
the present paper. 

For the diagram in Fig. la, four independent zero-threshold configurations exist, 
namely 

Case 1: one zero-2PT (e.g. m2 = rrir, = 0); 

Case 2: two zero-2PT's (mi = m2 = m4 = = 0); 

Case 3: one zero-3PT (e.g. m2 = m^, = = 0); 

Case 4: one zero-2PT and one zero-3PT (e.g. 1712 = = = = 0). 




Case 1 



Case 2 



Case 3 



Case 4 



Figure 2: Different zero-threshold configurations for the diagram in Fig. la. 

These cases are shown in Fig. |^ where bold lines are massive, narrow lines are massless, 
and the cuts corresponding to zero-thresholds are denoted by dashed lines. Note that, in 
case 1, having one more vanishing mass (e.g. or m^) does not produce a new zero- 
threshold configuration, and the corresponding cases can be considered together with case 
1, e.g.: 

Case la = Case 1 with ms = 0; 

Case lb = Case 1 with one more mass (not ms) equal to zero. 
The cases 1,2,3,4 are independent, and we shall discuss them separately. 
As in the paper we shall use general results on asymptotic expansions of Feynman 
diagrams. This subject was developed in several papers ||29| (see also ||3l| for a review). 



Before presenting these results in our case, we would like to discuss some issues con- 
nected with this approach. Let us first introduce the Taylor expansion operator 7^. Its 



^2PT and 3PT mean two- and three-particle thresholds, respectively. 
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action on the denominator can be represented as 



1 



E 



oo 



{v), {2{kp) - ey 



(4) 



[{k-py 



— m- 



j=0 



where 




r(z/ + j) 



(5) 



is the Pochhammer symbol. It is understood that this operator T acts on the integrands 
before the loop integrations are performed, and also that all the terms with equal powers 
of k on the r.h.s. should be collected together (e.g. the terms with k'^ and {kpY should 
be considered together). Below we shall also need a generalization of formula (D to the 
expansion in two or more momenta {k and one or two loop momenta). This generalization 
is straightforward: we can expand in all these momenta at the same time. We only have to 
remember that in this case the expansion is in the total power of these momenta, because 
inserting loop momenta in the numerator of the resulting massless integral effectively 
produces higher powers of k. 

Since the limit fc^ — > is not regular for the cases being considered, the true expansion 
cannot be obtained simply by applying eq. (D to all denominators where the external mo- 
mentum k occurs. If we did this, we would obtain infrared singularities in the coefficients 
of the expansion, even if the initial diagram were finite. This is due to the fact that in all 
zero-threshold cases we always need to let the momentum k go through at least one of the 
massless lines. Nevertheless, it is possible to correct this procedure by adding some extra 
terms. The general theorem on asymptotic expansions tells us how to construct this true 
expansion. Namely, the result is given by the naive Taylor expansion in small parameters 
(in our case, T^) plus terms given by Taylor expansions in a specific class of subgraphs (see 
below) which may also contain expansions in some of the loop momenta. So, this naive 
Taylor series produces infrared singularities while the extra series also involve ultraviolet 
ones, due to the appearance of loop momenta in the numerator. 

If the original diagram is finite, all the singularities should cancel^. Moreover, infrared 
and ultraviolet singularities should cancel independently. It is possible to write general 
formulae of asymptotic expansions in such an explicitly finite form, but such a procedure 
is more complicated. However, if we use dimensional regularization for both types of 
singularities, some contributions correspond to massless tadpoles and vanish. In this case, 
the formulae are much simpler, and we are going to use this approach below. It might seem 
in the resulting formulae that now infrared poles in e are cancelled by ultraviolet ones. 
The explanation is that omitting massless tadpoles effectively leads to a "re-distribution" 
of both types of singularities. So, from a practical point of view it is more convenient to 
apply the formulae where, formally, poles in e corresponding to infrared and ultraviolet 
singularities of the integrand are mutually cancelled. Moreover, cancellation of these poles 
for convergent diagrams will be a good check of our calculational procedure. 

Let us introduce the following notation: F is the original graph corresponding in our 
case to Fig. la, subgraphs of F are denoted as 7, and the corresponding "reduced graph" 
F/7 is obtained from F by shrinking the subgraph 7 to a point. Furthermore, is 

"^If it contains singularities of a specific type, those will survive. 



6 



the dimensionally-regularized Feynman integral with the denominators corresponding to 
a graph 7. For example, Jr corresponds to the integral (|l[). For our case, the general 
theorem yields 



r/7 



(6) 



where the sum goes over all the subgraphs 7 which (i) contain all the lines with large 
masses, and (ii) are one-particle irreducible with respect to the light (in our case, massless) 
lines. Tfc^ is an operator expanding the integrand in k and the loop momenta qi 
which are "external" for a given 7. The symbol "o" means that the polynomial in g^, 
which appears as a result of applying T to J-^, should be inserted in the numerator of the 
integrand of Jr/7- Note that although eq. (||) looks similar to the corresponding result 
for the large momentum expansion, eq. (2) of [^, there is an essential difference in the 
choice of subgraphs 7 contributing to the expansion. 




Figure 3: The subgraphs 7 to be included in the sum (H). 

Now, let us consider the contributions to the sum for the different threshold config- 
urations. We shall use the numbering given in Fig. la to indicate which lines are included 
in the subgraph: for example, {134} denotes the subgraph containing only the "left" tri- 
angle (lines 1,3,4) of Fig. la, {12345} = F, etc. These subgraphs are illustrated by Fig. ||, 
where bold lines correspond to massive propagators, narrow lines to massless ones, and 
dotted lines to the lines omitted in the subgraph (as compared with the "full" graph F). 
These dotted lines correspond to the reduced graphs r/7. We do not list the vanishing 
contributions containing massless tadpoles. 

Case 1. One 2PT is zero: m2 = = 0. Two subgraphs contribute to the asymptotic 
expansion: 
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(i) 7 = r 



(ii)7 = {134} 

1 1 

jTi in /T- /q\ 

^ ^[(A;-g)2]-2[g2].5 '^'^ [{k - pY - ml]''^ [{p - qf - mlY'^ [p^ - mlY^ ' ^ ' 

In cases la and lb we have the same contributing subgraphs (since the threshold config- 
uration is the same), and the only thing to change is to set = (case la) or 777,4 = 
(case lb) in the formulae (0)^(11) • 

Case 2. Two 2PT's are zero: mi = 1712 = = = 0, 1713 = m. In this case, four 
subgraphs contribute: 

(i) 7 = r ^ 

jn in /T-" . /q\ 

^[(p-g)2 -m2]-3[p2].4[g2].5 ^p_p)2].ip_^)2].2' y ) 

(ii) 7 = {134} 

^ ^ [(A; - qyp [g2]-5 [(A; - p)2]-i [(^ _ g)2 _ ^2].3 [p2].4 ' 

(iii) 7 = {235} contribution can be obtained from the previous one by the permutation 
1 ^ 2,4 ^ 5; 

(iv) 7 = {3} 

1 1 

jn in r-r- f^^\ 

^[(A;-p)2]-i[(A;-g)2]-2[p2].4[g2].5 [(^ _ g)2 _ ^2].3 " ^ ^ 

Case 3. One 3PT is zero: 7712 = ms = 7774 = 0. Two subgraphs contribute: 

(i) 7 = r 

(Tn era - Ti, - • (12) 

[{p - qyY^[p^Y^[q^ - mlY^o ^ [(^k - p^ - mlY^[{k - q^^^' ^ ' 

(ii) 7 = {15} ^ 

^ ^ [(A; - g)2]-2 [(p _ q)2Y-s [p2].4 ^^-p-^ p _ pf _ mf]-! [g2 - mlY^ ' ^^^^ 

Case One 2PT and one 3PT are zero: 7772 = fn-i = 7774 = 7775 = 0, ?77i = ?77. For 
this case, three subgraphs contribute: 

(i)7 = r ^ 

d"nd"o ^ Ti ^ • fl4) 

^^[(p-g)2]-3[p2].4[g2]z.5 ^P_p)2_^2].ip_g)2].2' ^ ^ 
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(ii)7 = {134} 



(iii)7 = {l} ^ 

Examination of eqs. (^-(|I6|) shows that, after partial fractioning is performed wher- 
ever necessary, the following types of contributions can occur in the expressions for the 
coefficients of the zero-threshold expansion: 

(a) two-loop vacuum diagrams with two (or one) massive and one (or two) massless lines; 

(b) products of a one-loop massless diagram with external momentum k and a one-loop 
massive tadpole; 

(c) two-loop massless diagrams with one (or two) powers of the denominators being non- 
positive. 

The algebra of contributions of type (a) is discussed in Appendix A, which is based on 
a slightly modified algorithm of |2^. The problem of numerators of massive integrals 



occurring in types (a) and (b) has been studied in p8| (see Appendix B of that paper), 
and we use the same formulae here. For the numerators occurring in the integrals of type 



(c), we have used a procedure based on a formula in and we explain the details in 
Appendix B. 

All the formulae (^-(|^), along with eq. (|^), can be used for any values of the space- 
time dimension n and the powers of propagators Vi. In the next section, we shall consider 
the results produced by the expansion (0) for some special cases. 



3. Analytic results 

To get analytic results for the coefficients of the zero-threshold expansion, we have 
written a computer program using the REDUCE system for analytical calculations 



Like the algorithm itself, it works for any (integer) powers of the propagators, and can be 
applied for both convergent and divergent diagrams. 

Many interesting physical applications require the calculation of the diagram in Fig. la 
with unit powers of propagators, the so-called "master" diagram. If all five Vi are equal to 
one (and fc^ 7^ 0) the corresponding diagram is finite as n ^ 4, and we shall calculate the 
corresponding results at n = 4. The algorithm also makes it possible to consider higher 
terms of the expansion in e = (4 — n)/2, or even results for arbitrary n. The individual 
contributions to the expansion ([7|)-(p!6D contain poles in e, which, for finite diagrams such 
as this "master" diagram, cancel in the sum (|^). This is one of the non-trivial checks on 
the algorithm. 

Let us write the corresponding "master" integral as 

00 

J(mi, m2, ma, m^, m^; k) = J(l, 1, 1, 1, 1; mi, m2, ms, m^, m^; k) = -vr^ ^ Cj {k'^Y , (17) 

j=0 
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j 


mOmmO 
(case 1) 


mOOmO 
(case la) 


mOmOO 
(case lb) 


OOmOO 

(case 2) 


mOOOm 
(case 3) 


mOOOO 

(case 4) 





2 

3 
2 


L 

-2 


L 

+C2-4 


-L^ + 2L 
-2C2 


-2C2 


-C2-3 


1 


24"-^ 
1 

16 




4 -'^ 

11/* 13 
+ 21,2 — 12 


+C2-f 


2 

-C2-i 


4 ^ n 2 

1 7 
— 2S2 — 8 


2 


1 r 

180 -'^ 

1 
180 


1 r 

90 -'^ 

_|_ 23 

' 540 


9 

_|_ 1 /• 673 
"T 3 ^2 1080 


3 9 

-IC2 + 1 


3 

2 A _L 13 
3^2 + 36 


5 T 9 

1 A 11 
"3^2 ~ 18 


3 


1 T 

1120 -'^ 

4^0 


1 r 

560 -"^ 

J ir_ 

1120 


1 r 

16 -'^ 

1 1 /• 1487 
"T 4 S2 3360 


/• 115. 

+ 2^,2 ^ 144 


5 T 

24 ^ 

— ^C^ -1- 763 
2 ^2 1440 


1 7- 2 _|_ 17- 

_!/■„ _ 2S7 
4 S2 576 


4 


1 r 

DoOO 

1 


1 r 

3150 
1 367 


1 !/■ 86939 


5 ^2 36 


2 r 1 13591 


— T77 / ; L 

_!/■„ _ 3841 
5 S2 9000 


10500 


1 63000 


5 ^2 252000 


5 S52 25200 


5 


1 r 

33264 J-^ 
1 


1 r 

16632 -'^ 
1 2531 


1 1 /■ 2828267 


17-2 17- 
18 

_|_ 1 A 973 

~r 3 S2 1800 


45 

1 /• 1 75683 


i r2 1 j_ 7- 

12 -t^ T 18 -t^ 
1 /■ 4049 


66528 


1 997920 


~r 6 S2 9979200 


3 S2 ~r 151200 


6 S2 10800 


6 


1 r 

168168 -'^ 
1 


1 r 

84084 -1^ 
1 4903 


1 1 /• 85046849 


— + ^-Z^ 

2 /• 1 139 
7 ^2 T 300 


13 7- 
210 ^ 

2 /• 1 2629903 


1 7'2 1 27- 

1 A 206741 


392392 


1 3923920 


"T 7 ^2 353152800 


7^2 1 5821200 


7 ^,2 617400 




compared 
with |17|, 
eq. (23) 


compared 
with H, 
eq. (105) 


L-terms 
checked 


compared 
with [|1|, 
eq. (104) 


L-terms 
checked 


compared 
with 0, 
eqs. (13)-(14) 



Table 1: Coefficients of zero-threshold expansion for the cases with only one mass. 



where, for zero-threshold configurations, the coefficients Cj can depend on logarithms of 

C, = \n\-k^) + Cf^ H-k^) + Cf\ (18) 

Powers of ln(— /c^) higher than two cannot occur in the cases considered, because the 
separate terms contributing to the expansion can have at most double poles in e. 

For cases when we have only one massive parameter m (all non-zero masses are equal), 
a number of exact results are known [|l^, 0, |T^; for some of them, closed expressions 



for the coefficients of the expansion are also available. They can be represented in terms 
of polylogarithms Li^v with < 3. While applying the algorithm described in Section 2 
we have considered all possibilities to have only one mass parameter in zero-threshold 
configurations. They correspond to the diagrams in the ffist column of Fig. 3 with all non- 
zero masses equal. The results for the ffist few coefficients of the expansion are presented 
in Table 1, where the first line gives the mass arguments of J, and the dimensionless 
coefficients (and other notations) are defined by 



J • i-'"(-;^j. c,-c(2) = -. (19) 

So, in Table 1 the following coefficients Cj (j = 0, . . . , 6) are presented, 

^ [m^y+^C, = cfL^ + cfL + cf , (20) 

and they agree with exact results wherever available (see the last line of the table). For 
cases lb and 3 exact results were not available, but we checked the logarithmic terms by 
use of the dispersion relations technique. 
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In cases 2 and 4, we have only one massive line (one mass parameter), and the expan- 
sion ([l^) is sufficient. However, for other zero-threshold configurations we can have two 
(cases la, lb, 3) or even three (general case 1) different masses. 

Let us start by considering case 1 with three different masses. We have obtained the 
following results for the lowest coefficients of the expansion (|1^: 



4(mf — ml) (ml — m3)(mf — ml) 



Co 

X I2 (in \-^] + In \-^] - 4 ] (mKml - ml) In ^ - ml{ml - mj) In ^ 



ml 

mt / \ mt I J \ ^ ~ ' Tn% — - ^2 



mi ( , m^ , m^\ „ o , m? , m? 



/ 2 2\/ 2 2\i lit ''''1 , ''4 \ 2/ 2 2\i I1 4 

+ (mi — mn)[m. — mo) In — ^ In — ^ + In — ^ — 2m J m, — mA In — ^ In — ^ 

^ ^ d/v 4 6J ^2 \^ 77^2 i 4/ ^2 ^2 

+2(m^ + ml){ml - ml)n{ml, ml) - 2{m\ + m^)(m^ - m^)7^(m|, m^) j , (21) 



(mf — mD'^lm: 



4 — m3)^(m 



m 



2^3 



X |2(m2-m2)(m2-m2) fin 




/ 2/2 2\2/ 2/ 2 I 2\ <-> 2 2\ i 

X m^ [m^ — mn) (m^ [m^ + m^) — im^mn) m — 5- 
V "^3 



2/ 2 

—m^yrrii 



m 



2\2, 



1712 



m^{m^ + ^4) — 2m-^^m^) In — 2 



-|-(m^ — m3)(m| — m3)(m^ — m|)(m3(m^ + m|) — 2m\m^^ 



I / 2 2\2/' 

+ (m^ — m3) (m; 



m 



2n2 



;m? + m|) fin 4 + In ^1 
\ mi mi) 



{m\ — m3)(m| — ml) In — | + 2ml{rn\ 

m,^ 



-ml) 



mt 



mi 



mS 



mi 



-|-2m3(m^ — m3)(m4 — m3)(m^ — m4)'^ In In 
-|-4mi(m^ 



m. 



((m?-mi)(m|-mi) 



m. 



2\2 



/ 2/ 2 2\2i 2/ 2 2\2i 

X m^[m^ — mr^) m — ^ + m^[mi — m^) m 



mq 



+2(mj^ — m3)(m4 — m^)(m^ — m4) (m-^m3 + m^m^ + mj^m4 — Smg) 
+2(m^ — mi)(m| — m^)^ (("^1 + 'm1){mi — m^) + 2m{ml{m\ — m^^^ H{ml, ml) 
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-2{ml - mlf{ml - ml) ((m? + ml){ml - ml) - 2m\m\{m\ - m^^ n{ml,ml) | .(22) 
In these formulae, 7i is a dimensionless function which is defined as 

H(™;,™=)^2L,(l-^)H-iwg). (23) 

It is easy to see that the function Ti. is antisymmetric in its arguments and therefore 
it vanishes when they are equal. This function is connected with the finite part of the 
two-loop vacuum integral with one massless and two massive lines (for details of this 
definition, see Appendix A). We have also obtained higher coefficients for this case (up to 
C3), but do not present them here because they are more cumbersome. 

If two of the masses are equal, mi = m^ = m, eqs. (|2l|)-(|22D can be simplified to 



•-^0 = ^^r^ ^ 2 m - mg - rrig In — In 

2 — m, r \ mi \ 



-ml In^ ^ - 4(m^ - ml) - 2ml ^("^^ ml) \ , (24) 



C, = ,,^.(J._^2)4 { (em^rn^ln^ + {m^ -ml){m^ -5m^ml-2mi)^ In 



9 

+3m^m3 In^ — 2m\{bm'^ — m\) In — 5- 
ml ml 



+ (m^ — m3)(m'* + ITm^mg + 2m3) + Gm^mg T-C{m'^, ml)^ . (25) 

In this case, we have got the coefficients up to C4. 

In another limit, m^ — » (this corresponds to the case la), we get for the lowest 
coefficients 

Co = / In 4 f-41 + In f-41 - 4 , (26) 
2[m{ — ml) ml I \ mfy \ mjl J 



- — „ „, „ ^ < mim^[mi + m^) In ^ + In ^ - 4 In ^ 

Amiml{mf-mi)^ 1 i 4v 1 4/ i i ^2 / I m|/ / m| 



-2m^m2(m? - ml) fin f— ^) + In ( -4]+ {mj - m^)4 , (27) 



mt / V mA 



?)Qm\m\{m'{ — mf)^ 



X l^m\m\{m\ + Amfml + ml) |^ln (^~^ j + (^"^ j ~ ^ j rr? 
-18mXK - <) (in + In - 4) 

4 4/ 2 2\2i 1 / 2 2\S/2 2\/4 22 4\l /\ 

— 2mim4(mi — m^) In — ^ + (m, — m^) (m, + m^) (m, — mim4 + m^) > . (28) 

mi I 



2 
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For this case (with two different masses), the exact result was also available in [|18|, and 
we have successfully compared our coefficients with this exact expression. 

Now, let us present the results for the coefficients (18) for case 3 (m2 = = = 0), 
when mi and are different: 



Co 



1 



-{ml + ml) 



In^^ + 

, mi 



2n' 



[m, 



m^ 



Himlml 



(29) 



^mfmKmf — m\) 



2m\ml{m\ — ml) 




— {m\ — ml){m\ + mg) In^ 



mi 



+ 



27r2 



- 2(m^ 



In ^ - 3 



m^ 



+ 2m\ml{m\ + mg) In 



,2\2 



(m2 + m2)7^(m2,m2)|, (30) 





?)Qmiml{m'l — m1)^ 



X ^mlml{ml-ml)\ml + ml) (hi [^-^ +^''(~^)) 

/ m"^ 9it'^ 

( 2 2 2\/ 4 2 2 4-\/i21 

— 3(mi — mg) (m^ + m^)\m^ — m^m^^ + m^) m — 2 H 



2 

0\/i 00/0 0\0 TXXi -] 

+3mim5 ( 3(mi — mg) +477117715(7711 — 7715) + 27r7 17775 1 m — 

' m^ 



^m^mv^im^ — m^im^ + 7^.5) [{m^ — 7^5) — 37n,i7n,5^ 

—Q{ml — ml)'^{m'l + mfml + m'^) H{ml^ml) >. (31) 



Again, we do not present higher coefficients because they are more cumbersome. For this 
case, we have calculated them up to C5. 

For all cases with different masses considered in this section (except eqs. (|26D -(^8D), the 
exact results were not available. Nevertheless, we have checked all the terms containing 
ln(— fc^) by use of dispersion relations. These logarithmic terms completely define the 
imaginary part of the results, which is non-zero for k'^ > 0. 

To see that the algorithm also works correctly for cases when the original integral is 
divergent, we considered an example corresponding to the diagram in Fig. lb, where the 
first and the fourth lines were taken to be massive (with equal masses) while the other 
three lines were massless. Effectively, this example corresponds to case 3 (one zero-3PT) 
with z/i = 2, 1^2 = t's = 1^4 = 1, 1^5 = 0. The exact result for this case was presented in |]TB[ 
(eq. (102)), and we have checked that our approach produces the correct coefficients for 
both the divergent and finite parts (we have calculated terms up to Cq). 



13 



Another possibility is to compare the obtained results with programs based on numer- 
ical integration. This will also illustrate the radius of convergence of the expansion. This 
comparison is considered in the next section. 

4. Numerical results 

In order to demonstrate the use of the zero-threshold expansion as a way to obtain 
approximate numerical results for self-energy diagrams, we shall discuss a number of 
mass configurations corresponding to diagrams that occur in the Standard Model. For 
this purpose, we take the masses of the W and Z bosons, and of the top and bottom 
quarks to be: 

Mw = 80 GeV, = 91 GeV, = 5 GeV, = 174 GeV. (32) 

The first example is a diagram containing a top-bottom loop to which two W bosons 
are attached, that contributes to the self-energy of the photon and the Z boson. The 
corresponding scalar integral is 

J{Mw,mh,mt,Mw,mh;k) . (33) 

Let us now neglect rrib (below, we shall see that this is reasonable when k"^ 3> ml), and 
consider instead 

J{Mw,0,mt,Mw,0;k). (34) 

It has one zero two-particle threshold and therefore it belongs to case 1. In Fig. ^ the 
approximations defined by 

Af 

j{N) ^ ^ ^^2y ^35) 

i=0 

are shown as curves, and, for comparison, values of J obtained by numerical integration 
[ P0| are shown as crosses. The position of the lowest non-zero threshold, at /c^ = 4M^ in 
this example, is marked by a vertical line. In the real part, we clearly see the logarithmic 
singularity as fc^ ^ 0. This is the region where, by construction, the asymptotic expansion 
works extremely well. For instance, at k'^ = 5000 GeV^, the N = 2 approximation is 
already just as good as the numerical integration, which is accurate to 4 digits, for both 
the real and imaginary parts. However, the expansion converges all the way up to the 
first non-zero threshold at k'^ = 25600 GeV^. At that point, the real parts of J^^-* and of 
J differ by about 15%, and the imaginary parts by 25%. 

We also examined the situation for k"^ < 0, where there are no physical thresholds and 
the imaginary part is zero. At small /c^, the behaviour of the real part for negative fc^ 
is quite similar to that for positive fc^, because it is defined by the ln|/c^| term in both 
cases. Moreover, because the signs of the terms alternate for fc^ < 0, the convergence of 
the asymptotic expansion is even a little better than for k"^ > 0. 

In Fig. ^, values of the integrals (|33D (dashed lines) and (Q) (sohd lines) are plotted 
to show the effect of neglecting rrify. At low /c^, their behaviour is obviously very different, 
the latter having a singularity at fc^ = 0, whereas the former has one at k"^ = 4:ml, but as 
fc^ increases above 4m^, they soon come close together. 
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Figure 4: The real and imaginary parts of J{Mw, 0, rrit, Mw, 0; k). 
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Figure 5: The real and imaginary parts of J{Mw, rrib, rrit, Mw, mb, k) and 
J{Mw,Q,mt,Mw,Q;k). 
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For momenta close to the Z mass shell, ~ 8300 GeV^, the error we make by 
neglecting m^, is of the order of 1%. The truncation error of the asymptotic expansion is 
already less than that if we use J'-^-'. 

The second diagram we studied contains a top-bottom loop with a W boson going 
across the inside. Neglecting mf, as before, it corresponds to 

J(mt,0,MH.,mt,0;A;). (36) 

This mass configuration also belongs to case 1. Below the first non-zero threshold, the 
behaviour of (|36|) is qualitatively very similar to the previous example, but because this 
threshold is now a three-particle threshold, it is different above the threshold. In this case 
the truncation error of J*-^-* on the threshold is only about 2% in the real part and 0.4% 
in the imaginary part. 

Finally, we study an example of case 3, 

J(M^, 0, 0, 0, Mz; fc) , (37) 

which corresponds to diagrams that contribute to, e.g., the electron or muon self-energy, 
if me or is neglected with respect to /c^. The results are shown in Fig. |^. The position 
of the first non-zero threshold, at /c^ = M^, is again indicated by a vertical line. For this 
case, we also checked that the numerical values (shown as crosses) agree with the results 
of other numerical programs [^, 

Contrary to case 1, the first term of the expansion, Cq, given by eq. (0), does not 
depend on ln(— fc^), and therefore, the integral has a finite limit as fc^ 0. At the 
threshold, where the truncation errors are largest, the errors in the real and imaginary 
parts of J^^^ are approximately 2% and 16%, respectively. 

In all the examples, the zero-threshold expansion provides approximations that are at 
least as accurate as numerical integration in a large part of the region of convergence, and 
can be evaluated much faster. 



5. Conclusions 

In this paper we have studied all possibilities for a two-loop self-energy diagram to have 
its lowest physical threshold at zero external momentum squared. Due to specific infrared 
singularities at /c^ = 0, constructing the small momentum expansion is more complicated 
than when all the thresholds are non-vanishing. Namely, instead of a regular Taylor 
expansion, we get a cut along the positive A;^ axis, starting with a branchpoint at A;^ = 0. 
As a consequence, our results in four dimensions contain ln(— /c^) terms, or even ln^(— A;^) 
terms if two physical thresholds vanish. So, technically this zero-threshold expansion is 
closer to the large momentum expansion than to a regular small momentum expansion 

For special cases, we have compared our results with known analytic ones. For some 
more complicated cases corresponding to diagrams with different masses occurring in 
the Standard Model, we made a comparison with the results of a numerical integration 



program based on the algorithm (for one case, we have also compared with ||24|| ). 



It shows that the zero-threshold expansion converges up to the first non-zero threshold. 
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Figure 6: The real and imaginary parts of J{Mw, 0, 0, 0, Mz, k). 
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and that, unless k"^ is very close to the threshold, only a few terms are needed to obtain 
accurate results. This comparison can also be considered as a check of the numerical 
programs. 

The algorithm constructed in the present paper covers the only "hole" remaining in 
the small momentum expansion of two-loop self-energy diagrams. Together with [p6| , 
it solves this problem completely. The large momentum expansion of these diagrams is 
described in I^S], and there are no remaining "holes" there. It is interesting to note that 



all these algorithms can be applied to two-loop three-point functions as well. 

The general problem of describing the threshold behaviour is still waiting for a solution. 
Here we have considered only a part of this problem, corresponding to the zero-threshold 
cases. In principle, the general theory of asymptotic expansions can be applied in some 
regions between the thresholds, whilst the description of the behaviour near the non-zero 
thresholds requires other approaches. 

A. D. would like to thank the Instituut-Lorentz, University of Leiden, and J.B. T. 
the Department of Physics, University of Bergen, for hospitality during visits when parts 
of this work were done. We are grateful to A. Czarnecki and J. Fujimoto for their help 
in comparisons with results based on other numerical approaches. A. D.'s research was 
supported by the Research Council of Norway. 



Appendix A. Two-loop massive vacuum diagrams 



To evaluate two-loop vacuum massive diagrams occurring in the separate terms (|^- 
(|16D contributing to the zero-threshold expansion, we mainly use the algorithms and 



formulae presented in [26, 28 



The integrals we are interested in are defined as 



/(i/i,i/2,z/3;mi,m2,m3) 



d> d"g 



1^2 



[{p - 



mill's ' 



(3^ 



The general case of this integral with arbitrary n and z/'s was considered in E^. By use 



of the technique a result in terms of hypergeometric functions of two variables was 
obtained. For the case when all z/'s are equal to one and n is arbitrary, this result was 
reduced to hypergeometric functions 2F1 of one variable (see also in [l^). Expanding in e, 
we get standard results in terms of dilogarithms or Clausen's function given by eqs. (4.9), 
(4.10) and (4.15) of ^ (see also in |ll], H). 

In this paper, however, we always have at least one of the masses in (|3^) equal to zero. 
Generally, this simplifies the calculations but requires some changes in the algorithms as 



compared with the general mass case 



If one of the masses (e.g. ms) is zero, the result for the integral 
in terms of a Gauss hypergeometric function. 



can be expressed 



J(z/i,z/2,i^3;"^i,"^2,0) = 7r"z '"(-"^2) 



2\n~i/i—i/2—iy3 



^^3 



r(z/i + 



^3 



r Z/2 + Z^3-^ r(z/i + Z/2 + Z/3 



n) 



T (z/i) F (z/2) F (f ) F (z/i + + 2z/3 - n) 
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X 2^1 



Z^l + - f , + ^2 + 1^2. - n 

+ + 2z^3 - n 



mil 



(39) 



where it is understood that i (— m = (m This result can easily be obtained using 



the formulae of ^4], The symmetry (mi,z/i) ^ {^121^2) can be seen by use of the 
well-known transformation 



2F, \' 



z - 1 



(40) 



For the cases when mi = 1712 or m2 = the formula (^) gives the correct answers: 
/(z/i, z/2, z/3; m, m, 0) = 7r"i2-2n(_^2)n-^i-'^2-'^3 

r (f - z/3) r (z/1 + z/3 - 1) r (z/2 + ^3 - f ) r {u, + u2 + iy3-n)_ 



X- 



r(z/i)r(z/2)r(f)r(i.i + i/2 + 2z/3- 



(41) 



/(l/l, 1/2, i^3; m, 0, 0) = 7r^2-2n(_^2)n-.:-.2-.3 

r (f - Z/2) r (f - z/3) r (1/2 + ^^3 - f) r {u, + u2 + u,-n) 



r(z/i)r(z/2)r(z/3)r t 



These results are well-known, and we present them here for completeness only. 

For the case when all z/'s are equal to one and n — s> 4 (e — 0), the result for 
be written as 



/(1, 1, l;mi,m2,0) = vr 



4-2e 



2(i-e)(i-25) 



X < — yif^i + f^l) "I — ("^1 In(m^) + m2 ln(m2)) 



e 



-2{ml In^(mi) + ml ln^(m2)) + -{ml + ml) In^ ^ + (m^ - ml)7i{ml, ml) 

2 777- 

where 7i is a dimensionless antisymmetric function defined as 

HK,»^,^2Lji-;|)4wg). 



(42) 



can 



(43) 



(44) 



It is connected with the function F used in (see Appendix B, eqs. (B.3)-(B.7) of |2^ ) 
via 

^ F(m'i,mn,mVi 1, m? / m? , 
V 1) 2, In^ + ln-^ 



H{mi,m2 



lim 

ma— *0 



2 2 
m| — 7712 



2 ml I ml 



(45) 



We cannot use the original function F here because it has a logarithmic singularity as 
m3 0. 
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We need to calculate integrals (^) with arbitrary (integer) powers of denominators. 
For the cases when we have one mass parameter only, it is convenient to use the simple 
explicit formulae (Jl), (|2). If we have different masses, however, it is better to use 
recurrence relations obtained by using the integration-by-parts technique (in a way 
analogous to |0). For the case of the integrals (0), these relations were considered in 
26]. The solution to these relations presented there (see eq. (5.3) of p6|) contains masses 



in the denominator. This can be avoided by employing the connection (5.4) of the same 
paper. Note that for = the determinant of the system of equations corresponding 
to the recurrence relations becomes A(m^,m2,0) = —{mf — rn^"^. 

Of course, application of these recurrence relations is equivalent to using connections 
between contiguous hypergeometric functions 2-^1 (see, e.g., in ||3^) corresponding to 
eq. (|39|). Furthermore, in cases when we have close values of the masses, it may be 
better to use the explicit result (|39D to expand in the mass difference up to the order 
corresponding to the required accuracy. 



Appendix B. Two-loop massless integrals with numerators 



Here we shall consider massless integrals 



^^2, ^^3, ^A. y^) = J ({^^^}; 0, 0, 0, 0, 0; k) 



(46) 



where J corresponds to the original integral (|l]). The main properties of these integrals 
are well-known, and many of them can be obtained by use of the integration-by-parts 
technique (we had collected some of them in Appendix A of [^). In particular, 
when any of u's is equal to zero, the corresponding integral can be calculated in terms of 
gamma functions, for any value of the space-time dimension n. If all i/'s are equal to one 
and n = 4, the result is well-known (— 6C(3)7r^/fc^). For higher powers of propagators, 
integration by parts [^] can be used, which was implemented in programs like [39 . 



In this appendix, we would like to discuss a "non-standard" way of calculating the 
integrals in cases when some of the z/'s are negative. So, in fact, we get the corresponding 
denominators in the numerator. This approach is based on the formula presented in 
ref. [0 and was used in this paper. The general case of this formula can be applied to 
any one-loop integral with tensor structure in the numerator. As a result, one gets a sum 
of tensor structures multiplied by scalar integrals in other dimensions. We shall need only 
a simple corollary of this general formula below. 

In fact, we have two different cases to consider: (i) 1/3 < and (ii) any other u (for 
example, u^) is negative. Let us consider the first possibility, the second one can be treated 
analogously. If z/3 < 0, this corresponds to getting 



P3 



p"^ + - 2{pq) 



(47) 



in the numerator. After expanding and cancelling p"^ and against the corresponding 
denominators, we get integrals of the form 



J J d> (Tq 



(pq) 



M 



[pT [qT {{k - p)T [{k - qy] 



(48) 
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Now, let us write (pq)^ as p'^^ . . . p^^' ■ ■ ■ %m and consider one of the subloops (an 
analogy with 0] can be noticed). According to ||32|, it can be written as | 

X E (^)'{b]1*=]'=}^,.,,^„rg-..+A+.)rg-..+A)r(..+.,-^-A).(49) 

A, K 
2X+K=M 

where \ [g]'^[k]'^\ denotes a symmetric tensor structure each term of which is con- 

structed from A metric tensors gf_i,fij and k vectors k^., the coefficient of each term being 
equal to one. It is easy to show that such a tensor structure contains 

2^ A! k\ ^^^^ 

terms. Moreover, when contracted with p^^...p^'^', each term gives the same result, 
namely {p"^)^ [kp)^. 

So, now we obtain a sum of the terms of the form 



d> ^ , (51) 

[pT [{k-pyr 



with known coefficients (^). Finally, the scalar products [kp) in the numerator can be 
represented in terms of the denominators and k"^, and the resulting integrals are trivial. 
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